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The folded node is a singularity associated with loss of normal hyperbolicity in systems 
where mixtures of slow and fast timescales arise due to singular perturbations. Canards are 
special solutions that reveal a counteractive feature of the local dynamics, namely trajectories 
that flow from attractive regions of space into repulsive regions. An alternative way to 
model switches between timescales is using piecewise-smooth differential equations. There is 
presently no adequate theory or method for relating slow-fast and piecewise-smooth models. 
Here we derive the analogous piecewise-smooth system for the folded node by pinching phase 
space to sharpen the switch between timescales. The corresponding piecewise-smooth system 
contains a so-called two-fold singularity, and exhibits the same topology and number of 
canards as the slow-fast system. Thus pinching provides a piecewise-smooth approximation 
to a slow-fast system. 


The purpose of this paper is to study the relation between discontinuities and singular per¬ 
turbations in dynamical systems, by focusing on their singularities. We do this by forming a 
piecewise-smooth model of an important singular perturbation problem using the method of pinch¬ 
ing introduced in j^. Pinching replaces a slow-fast timescale separation with a discontinuous switch, 
replacing Fenichel’s slow manifolds jl^ with Filippov’s sliding surfaces [l7|. In particular, it pro¬ 
vides a piecewise-smooth model that captures, both qualitatively and quantitatively, the intricate 
dynamics of the smooth singular perturbation problem. The study of such singularities in discon¬ 
tinuous systems predates those in singularly perturbed systems % Il6l . |24| , while the latter has seen 
more progress in the study of oscillatory dynamics, see e.g. js]. Here we show that the analogy 
between singularities in the two types of system is more than superficial, paving the way for a more 
rigorous study of piecewise-smooth systems as a means to studying singular perturbations in the 
future. 
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I. INTRODUCTION 


Highly nonlinear changes in the dynamics of a system can be modelled by differential equations 
whose solutions vary rapidly near certain thresholds. Whether those jumps are smooth but fast 
varying, or are truly discontinuous, changes the way they are treated mathematically. It is natural 
to assume that a limit exists in which ‘fast-but-smooth’ becomes truly ‘discontinuous’, but how to 
characterise that limit is by no means obvious. 

Analytically this might not be surprising, because in many cases the limit may be singular 
(for an interesting aside on singular limits see ^), meaning that solutions of the smooth system 
do not limit to the piecewise-smooth system in a regular way. Computationally the difficulty 
in studying the nonsmooth limit of smooth systems lies in the fact that, by their very nature, 
they become extremely stiff and numerical methods fail to converge in the limit of interest (l^ . 
As a result, the relation between smooth and piecewise-smooth dynamical systems theory remains 
poorly understood. The recent growth of piecewise-smooth dynamical systems theory has produced 
numerous forms of discontinuity-induced singularities, bifurcations, and chaos. It is interesting to 
ask whether these have any counterpart in smooth systems, where novel behaviours have also been 
attributed to abrupt change. For this purpose a relation between smooth and piecewise-smooth is 
necessary. (Examples of the discontinuity-induced phenomena of interest are sliding bifurcations 
|3, 3, explosions jl^, grazing singularities [l7, 27|, and non-deterministic chaos [71. Examples 
of the smooth system phenomena of interest are canard explosions [ 4 1 and mixed-mode oscillations 



], ideas from singular perturbations, nonsmooth dynamics, and nonstandard analysis, were 
combined to develop a method called pinching, which approximates a fast-but-smooth change by 
a discontinuity. Pinching characterises the dynamics in the nonsmooth limit at least qualitatively, 
and, as we show here, even quantitatively, by preserving certain singularities and associated geom¬ 
etry. The present paper investigates the method by applying it to a known system with slow-fast 
dynamics, characterized by an invariant manifold of slow dynamics which loses normal hyperbol- 
icity, taking the form of the so-called folded-node system in a singular limit. A relation between 
so-called eanard phenomena in smooth and piecewise-smooth systems is derived. 

Pinching is a way of deriving piecewise-smooth models that capture key geometry of singular 
perturbation problems. While the method so-far developed does not provide a rigorous approxi¬ 
mation in the analytic sense, it faithfully captures singularities and bifurcations that arise when 
a smooth system suffers rapid change, encapsulated in a piecewise-smooth differential equation. 
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capable of providing accurate estimates of bifurcation parameters, as shown previously in a study 


0 . 


of the van der Pol system in 

In the remainder of this section we introduce the canonical model used to study the folded node, 
a set of ordinary differential equations with a singular perturbation parameter e. In section ITII we 
apply pinching to obtain a piecewise-smooth system dependent on e 7 ^ 0 , whose phase portrait 
resembles the smooth system’s singular limit e = 0. Neither the pinched system, nor the singular 
limit of the smooth system, accurately represent the dynamics found when e is small but nonzero. 
In section m we motivate pinching better by first making an exponential rescaling of the phase 
space, yielding a system with similar qualitative features to section m we then give numerical 
evidence, by means of boundary-value problem continuation, of the continuum of canards emerging 
in the nonsmooth limit. In section nYi we show how this can be improved, by shifting the focus of 
the exponential rescaling, and pinching again. This second approximation possesses more intricate 
dynamics, which we show is in one-to-one correspondence with solutions of the smooth system for 
e 7 ^ 0. We make some closing remarks in section 0 including a simple smoothing of the pinched 
system from section [ITl that shows how the number of canards in a flow increases, through a series 
of bifurcations, as a change at some threshold tends towards discontinuity. 

We begin with a slow-fast system whose key features can be considered fundamental to the 
understanding of the canard phenomenon. The system is a set of ordinary differential equations 
for two slow variables, x and y, and a fast variable, 2 . (Some features of these systems can be 


generalised to arbitrarily many fast and slow variables, see 


29|). The timescale separation is 


introduced by a parameter £ satisfying 0 < e <C 1. A general such system can be written 


X = gi{x,y,z;£) , 

y = 92{x,y,z-,£) , (1) 

ei = h{x,y,z-,£) , 

in terms of smooth functions h,gi, and g 2 , where the dot denotes differentiation with respect to 
the (slow) time t. For \h\ > e, it is easy to see that i is much larger than x and y, so on the fast 
timescale solutions are attracted to, or repelled from, an e-neighbourhood of the null surface of the 
fast variable, 


:= {h{x,y,z,£) = 0} , (2) 

known as the critical manifold] see figure [H As shown by Fenichel [l^ , within that neighbour¬ 
hood he (generally non-unique) perturbations of the critical manifold, called the slow manifolds, 
hypersurfaces which are (locally) invariant in the flow, and are hyperbolically attracting or repelling 
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provided that dhjdz 7 ^ 0. The notion of local invariance means that solutions can leave a slow 


manifold only at its boundary, if one exists 


m. 



Figure 1: Critical manifold C°, and dynamics away from the fold (the j/-axis). Solutions evolve quasi-vertically on the 
fast timescale (double arrows) to/from the e-neighbourhood of the critical manifold, where the slow timescale dominates 
(single arrows). 


Normal hyperbolicity of the critical manifold is lost if dh/dz vanishes. This scenario occurs 
generically in three dimensions when there is a fold in with respect to the flow (see figure [1]) . A 
fold is a set of points where the conditions 


^ = = ( 3 ) 

are satisfied, with the inequalities ensuring that the critical points along h = 0 are indeed folds, 
and not higher order degeneracies. The flow’s projection onto the x-y plane is transverse to the 
fold except at points where 


d 0 

h = {91,92)- (Tr-,^)/i = 0 


The three conditions in equation 


dy 

and equation ([5]), subject to non-degeneracy conditions 


( 5 ) 


define an isolated point, which we can place at the origin (x, y, z) = ( 0 , 0 , 0 ) for e = 0 , where the 
system can be transformed into the local canonical form js, 26| 


X = by + cz + h.o.t. 


y = 1 + h.o.t. 
ez = X + z^ + h.o.t. 


(6) 
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This provides a local model for the system of interest, and from now on we omit higher order terms 
(h.o.t.). In doing so we must take care, for we can scale e out of the leading order terms by sending 
{x, y, z, t} {xe, yy/e, Zyfe^ (following | 3 |), or {x, y, z, b, c} {xe^, y, ze, 6 e^, ce}, and while 
the result is a formally correct leading order system, the e-orders of the new variables (x, y, z, t) must 
be considered in any subsequent approximation. Since the order of certain quantities is essential 
here, we do not make such a substitution at this stage. 

A. Projection onto the critical manifold 

This paper presents a method for characterising slow-fast dynamics for a system with small e, 
without taking e to zero. The more common method of studying a system such as jbl) involves first 
taking the singular limit, e = 0, and subsequently considering e > 0 as a perturbation. Certain 
features of the singular limit will appear in a different guise (for e 7 ^ 0 in fact) in our nonsmooth 
approach later, so let us review these hrst. 

Setting e = 0 (and neglecting higher order terms) reduces system (jb]) to a differential-algebraic 
equation 


X = by + cz , 


y = 1 , 

0 = X + z'^ , 

called the reduced system, see e.g. [s, 
which h = X + z‘^ =0. This flow therefore satisfies h = 0 which, using system jb]), gives 

by + cz 


2b|. Its solutions are restricted to the critical manifold C^, on 


0 = h = by + cz + 2zz 


z = -- 


2z 


Combining this with the second component y = 1 from system ([b]) yields a dynamical system on 
given by 






( 7 ) 


This system is undefined at y = z = 0 but, importantly, arbitrarily close to this point the righthand 
side of equation (j7]) is generally nonzero. The local phase portrait is determined by the 2 x 2 in 
equation jT]), whose trace and determinant are respectively 


A 1 A 2 = 2b and Ai -|- A 2 = c 


(8) 
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where Ai ^2 are solutions of the characteristic equation 0 = A^—cA + 26 assigned such that |Ai| > IA 2 I. 
Let n = A 2 /A 1 . The singularity is then classified (see for example 


26|) as: 


a folded node 

if 

^ > 0 , 

( A 1 A 2 > 0 ) 

(9) 

a folded saddle 

if 

y<0, 

( A 1 A 2 < 0 ) 

(10) 

a folded focus 

if 

fj, (Ei C . 

(Ai=A^ ) 

(11) 


The three cases are illustrated in hgureEl Their names reflect the fact that, if we omit the prefactor 
— l/2z from equation ([7]), the remaining linear system 


y 

z 





has an equilibrium at the origin, and this is a node if /r > 0, a saddle if /U < 0, or a focus if /r G C. The 
term ‘folded’ is required because this linear system is topologically, but not dynamically, equivalent 
to equation ([7]), being obtained from it by a time scaling 1 1 -^ This scaling changes sign with 

z, and is singular at z = 0. As a result, the dynamics is similar in z < 0 on the attracting branch 
of the critical manifold , and is similar up to time-reversal in z > 0 on the repelling branch of C^. 
The singularity of the scaling at z = 0 means that the equilibrium of the linear system is not an 
equilibrium of equation ([7]), instead referred to as a folded equilibrium^ which, unlike an equilibrium, 
solutions can cross through in finite time. Solutions of equation ([7]) that cross through y = z = 0 
are shown in figure O with inhnitely many in (i), only two in (ii), and none in (iii). Those that 
pass from z < 0 to z > 0 through the folded node or folded saddle correspond to canards, which 
are discussed in the next section. 



The classification into folded node/saddle/focus, for different values of /i, is also used to classify 
the system ([6]) for e nonzero. The remainder of this paper will be concerned solely with the 
folded node case, /r > 0. Substituting y = A 2 /A 1 and equation ([8]) into equation ([6]), then scaling 
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(x, y, z,t,e) ^ (x/Af, y/Af, —z/Ai, t/Af, —e/Xf), gives to leading order 


y = I, ( 12 ) 

ei = X + . 

Since much of the foregoing analysis will involve a natural length scale \h\ = £, for convenience we 
define a new function 


u{x,y,z-,e) = h{x,y,z]£)/£ , (13) 

then in the variables {u,y,z}, the system m becomes 

£u = ^y—{l + iJ,)z + 2zu, 

y = 1 , (14) 

z = u . 

The remainder of the paper is a study of this system for /r > 0 and e > 0. Before we apply pinching 
to this system, we must state some preliminaries concerning the so-called canard type solutions that 
make it so interesting. 


B. Canards 

Throughout this paper we define a canard as a solution that evolves from an attracting invariant 
manifold, to a repelling invariant manifold, via some singularity that faciliataes the transition. In 
the traditional setting of a smooth system with slow and fast timescales, the invariant manifolds are 
surfaces of slow dynamics, as we introduce in equation (1151) below. In the setting of discontinuous 
systems as we introduce in section El the invariant manifolds are regions where solutions slide 
along the switching manifold (the discontintuity set). Strictly speaking, we use the term ‘canard’ 
exclusively for maximal canards, which are the solutions of the types above that spend the maximum 
possible time on the repelling invariant manifold;in the setting described in this paper, this time 
can be regarded as infinite. 

For real y and noting |/r| < 1 by definition, the ratios zjy = y/2 and zjy = 1/2 are satisfied 
by the weak and strong eigendirections associated with the folded equilibrium of equation ([7]). The 
solutions along these directions are canards, known as the weak and strong singular canards of the 
singular (e = 0) system. 
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The non-singular (i.e. e 7^ 0) system (I14p also has two particular solutions satisfying z/y = /i/2 
and z/y = 1/2. We label these = {u{t), y{t), z{t)}, and solve equation (fT41l to find 


7”‘(t) = {f, t. f* } . = 


(15) 


These solutions are canards, since as t —>■ ±00 they lie in an e-neighbour hood of the attracting 
and repelling branches of the critical manifold, implying that they tend towards attracting and 
repelling slow manifolds. They form simple curves that are e-close to the weak and strong canards 
of the singular system, therefore is called the weak primary canard, and 7** the strong primary 
canard, (though where possible without ambiguity we omit the word ‘primary’). 

It is clear from the local phase portrait (see figure [2](i)) that, besides the primary canards, 
the singular system contains a whole family of canard solutions through the singularity, forming a 
continuum between the weak and strong solutions. For the non-singular (0 < e <C 1) system there 
may exist a number of other canards, termed secondary canards, though unlike the singula r sy stem 
they will generally be finite in number, and of a more complicated topology (see for example |lll.l28l|i. 
Secondary canards have been shown [h, ^ to rotate around the weak (primary) canard near the 
origin, and to asymptotically align with the strong (primary) canard as t —)• ±00, as sketched in 
figure O They are neither easy to express in closed form, nor easy to simulate numerically. To 
study secondary canards analytically, Wechselberger [28| applies a parameter blow-up, then moves 


to cylindrical coordinates centred on the weak canard, and takes the variational equation along 
the weak canard to obtain a Weber equation, whose solutions describe small oscillations that the 
secondary canards perform around the weak canard. In this paper pinch the e-neighbourhood of the 
critical manifold onto the manifold itself, thereby deriving a piecewise-smooth system that exhibits 
analogous behaviour. 

For e nonzero, a first approximation for the slow manifolds is that they lie in an e neighbourhood 
of the critical manifold (at least where it is normally hyperbolic). A better approximation, and 
one we will use later, is to note that the slow dynamics lies not on, but is stationary with respect 
to, the critical manifold, and hence lies close to (actually in an e^-neighbourhood of) the nullcline 
u = 0. Solving M = 0 in equation (1141) gives the surface 

/r + 1 


Vy = <{u,y,z) : u = 


Az 


(16) 


These two approximations (e-close to C® and e^-close to Vy) must be consistent, that is, the sur¬ 
face Vy can only approximate a slow manifold where it lies in the e-neighbourhood of the critical 
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Figure 3: (i) Primary and secondary canards in the neighbourhood of the critical manifold. Secondary canards connect 
an £-neighbourhood of the attracting and repelling branches of the critical manifold, their tails lying parallel to the strong 
(primary) canard 7 “*. Near the singularity they rotate around the weak (primary) canard 7 ™*. (ii) The flow circulating 
around the weak canard: simulation of equation (1141) in the plane y = —1. 

manifold. Near z ~= 0 this cannot hold, while for large z the surface Vy is approximated by 

^0 = I , (17) 

and for this to lie within e of the critical manifold, which corresponds to the region |u| < 1, we 
must clearly have |(^ + 2)/2| < 1. Combining this with the folded node condition /r > 0 gives the 
range of permitted values for the parameter as 

0 < ^ < 1. (18) 

This motivates the choice we made earlier of defining // as A 2 /A 1 rather than its reciprocal. As 
an illustration of the condition above, note that the weak and strong canards are two particular 
solutions that lie on the slow manifolds. These he on Vy because their distances u from the critical 
manifold are fixed at u = /r/2 and u = 1/2, and moreover |^| < 1 guarantees that both of these lie 
in the required neighbourhood, since |/u/2| < 1 and |1/2| < 1. 

II. PINCHING 

Pinching, at least in the form used here, was introduced in Q]. To illustrate the method, we 
first demonstrate it in its crudest form. Essentially, we assume the system is dominated by fast 
dynamics for |n| > 1, which we leave untouched. We assume that slow dynamics dominates in 
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|tt| < 1 , in such a way that dynamics in the u direction can be neglected, so we collapse the entire 
neighbourhood |rt| < 1 onto the critical manifold tt = 0 . 

The ‘pinch' is enacted by a piecewise-smooth transformation of the variables, in this case intro¬ 
ducing a new variable 

U = u — sign(tt) , (19) 

for |m| > 1, and omitting |u| < 1. Equations m then become a piecewise smooth system 

eil = - {fi + l)z+ 2z s\gn{U) + 0{zU) , 

y = 1, 

z = U + sign(C/) . 

where 

' U/\U\ if C/ 7 ^ 0 , 

[- 1 ,+ 1 ] if u = o. 


( 20 ) 


sign([/) G 


( 21 ) 


The dynamical theory of such differential inclusions was described by Filippov |l7l |. For U 0, 
equation (j20jl specifies the fast dynamics uniquely. On U = 0, called in piecewise-smooth dynamics 
the switching manifold, the righthand side is set-valued. The dynamics it gives rise to, however, is 
rather simple to describe. 

Consider a point p on the switching manifold, so U\p = 0. If [/ is nonzero there and its sign does 
not change with the sign of U, then p is the arrival point of a solution lying on one side of C/ = 0, 
and the departure point of a solution lying on the other side. Concatenating the two solutions gives 
a unique, continuous but non-differentiable, solution, that crosses the switching manifold at p. This 
therefore takes place on [/ = 0 where Az^ < — {pi+ 1 )^^)^ (the ‘bow-tie’ region in figure 0 ]). 

If U changes sign at p then the flow cannot cross through the switching manifold there. Two 
solutions of equation ( 1201 ) meet at p, each arriving from either side of 17 = 0 (the lower re¬ 
gion on U = 0 in figure 0|) , or each departing (the upper region on U = 0 in figure 0]) . They 
occupy the regions on t/ = 0 where — (/r-|- 1 ) 2 ;)^; the region where z < 0 and 

4z^ — {^y ~ (/^ + l)'^)^ > 0 attracts the flow outside the switching manifold, while the region 
where z > 0 and 4z^ — ~ (/^ + 1)'2)^ > 0 repels it. These are the pinched analogues of a 

point p on the attracting and repelling slow manifolds, respectively, of system (I14p . To find the 
flow on them we solve as we did for the slow flow projected onto u = 0 in section fixing 
0 = U = — (p + l)z -|- 2zz. Solving for z gives dynamics on 17 = 0 defined by 

^ ^ % -2 

(22) 

f -ih + l) 



l=—f 

u 

1 2z \ 



11 


We say that system (1221) defines sliding dynamics on [/ = 0, and its solutions are known as sliding 
orbits. This is clearly sensible as an analogue of the slow dynamics, since it is equal to equation 
which describes the original system projected onto the critical manifold. 

Finally, notice what happens at the boundaries of the regions of sliding and crossing, where 
U = 0 and — (/r + l)z ± 2z = 0. The field in [/ > 0 is tangent to the switching manifold along 

— (/X + \)z + 2z = and curves away from the manifold since it satisfies U = ^ — (/U + l) + 2 = 
1 — ^ >0. Thus the flow in [/ > 0 at such points carries a solution away from the switching 
manfiold. The field in 17 < 0 is tangent to the switching manifold along ^y — (y + \)z — 2z = 0, 
and curves towards the manifold since U = ^ + (/x + l) + 2 = 3(^ + 1) > 0, so solutions can only 
enter the switching manifold at such points. 

All this gives the simple dynamical portrait shown in figure |4l A feature of this system is that it 
contains solutions corresponding to the weak and strong canards, and a continuum of other canards 
between them. Evidently, recalling section di the sliding dynamics corresponds to represent the 
continuum of canards that exist in the smooth system’s singular limit e = 0. But the pinched 
system also allows crossing of the switching manifold, which the singular (e = 0 ) smooth system 
does not, so it can only apply to e 7 ^ 0. It will emerge that the correct canard structure can be 
captured by modifying the pinch slightly (looking ahead, we do this by pinching around the surface 
(I17p instead of the critical manifold). In section UTTl we re-consider the choice of coordinates above, 
inserting a preparatory step (a “microscope") that more fully motivates the process of pinching, 
before changing the focus of the microscope and the pinch to capture more precisely the secondary 
canard structure in section nYi 




Figure 4: Pinched system, in {U, y, z) coordinates (left) and corresponding (X, y, z) coordinates (right) if X = a; —sign(a;). 
The upper and lower regions are the repelling and attracting sliding regions, separated by a ‘bow-tie’ (shaded) made up 
of crossing regions. Double arrows indicate trajectories outside the surface, single arrows indicate sliding trajectories. 
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III. FIRST APPROXIMATION: A CONTINUUM OF CANARDS 

The microscope, introduced in the context of nonstandard analysis j^, is an exponential scaling 
of variables that attempts to resolve the interaction between slow and fast dynamics. We use it 
here to motivate the pinch that will follow. Let 

V = := |u|^sign u , (23) 

in terms of which equation (I14p becomes 

i) = u {2z + (Ly — (^ + 1 ) 2 ) , 

y = 1, (24) 

i = . 

This is simulated in figure [5Ki). The weak and strong canards now satisfy = zjy = /i/2 and 
^i/b] = zjy = 1/2 respectively (by direct calculation). The projection onto v = 0 gives the same 
y-z system (12211 (or equivalently ([7])) as before the microscope. 

1 

For |i;| < 1 and e <C 1, the flow of system (I24p is dominated by the term v FI in the equation 
for V, and so lies close to a set of fibres with constant y and z, connecting the surfaces v = ±1. To 
approximate this, we pinch the two surfaces u = ±1 together using the transformation 

V = V — sign V (25) 

for jul > 1 , giving the system 

V = (U + signU){ 2 z+(fy-(/i + l)z)(U + signU)-V[d} , 

y = 1, (26) 

i = (l/ + signU)V[^l , 

as simulated in figure [5Kii). Note that prior to pinching, both the weak and strong canards lie inside 
the region |u| < 1, therefore they are not part of the system (|26p for U 7 ^ 0. Instead, they now lie 
on the switching manifold U = 0. At the switching manifold, V = 0, equation (1261) reduces to 

V = ^y - {n +l)z + 2z s\gnV , 

y = 1, (27) 

i = signU , 

which is equivalent to system ()20p up to an £ scaling in y,z,t. The crossing and sliding dynamics 
on the switching manifold are therefore exactly as described in equation ([22]) for system (I20p . The 
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microscope, given by the transformation u>-^ v, provides the motivation for the pinch u^V. The 
result is a system of ‘fast’ dynamics for V ^ Q given by equation (I26jl . and slow dynamics that 
either crosses the sliding manifold, or slides along it as described in equation (j22p for system (j20p . 
An immediate consequence is that, similar to the previous section (where the pinch was applied 
without a microscope), pinching here gives a continuum of canards solutions. In the next section 
the microscope is used more powerfully, to resolve the different canards. 



Figure 5: The microscope and pinch, (i) Flow in the microscope system II24II simulated in the plane y = —1 with e = 0.05 
and fj, — 1/8.5. (ii) Flow in the system II26I) obtained by pinching together the surfaces v = ±1 in (i), with the vertical 
axis rescaled for clarity. The primary canards are indicated. 



Figure 6: Sliding flow after the pinch, given by equation (1221) in the regions > (^y — (p + l)z)^ on F = 0. Showing 
the strong (bold) and weak (dashed) primary canards. Orbits to the right of the strong canard in the upper half space 
(the repelling sliding region), and to the left of the strong canard in the lower half space (the attracting sliding region), 
form a continuum of canards. 
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A. Computations in the nonsmooth limit 


Before resolving the secondary canards by adjusting the pinch, one may immediately ask whether 
the original canard structure is restored when we smooth out the piecewise-smooth model (I20jl . The 
obvious way to achieve this is to replace the sign function with a sigmoid function. Neglecting the 
0 {zU) term in equation (I20p . we directly replace sign with 

sign(t/) !-)• tanh(fc[/) , 

for large positive constant k, and analyse the resulting system numerically using a technique based 
on the numerical continuation of parametrised families of two-point boundary-value problems as 
presented in [ll|. The piecewise-smooth system (1201) is obtained in the limit A: —)• oo, and the 
original smooth system da is regained as A: —)• 0. 

The computations below are all made for ^ = 1/8.5. Figure [7] shows attracting and repelling slow 
manifolds for different values of k, found by computing solutions that pass between lines chosen on 
the attracting and repelling branches of the critical manifold. Canards occur where the attracting 
and repelling branches intersect transversally. The magnified images show an increasing number of 
rotations around the weak canard, and an increasing number of intersections (i.e. canards), with 
increasing k, that is, as the smoothing function tanh(A:?7) approaches sign([/). 

Figure m shows ten branches of canards ( 7 ^ to 7 ® and 7 ®*) continued as the stiffness parameter 
k varies (see {h| for details on such continuation for canards). Each point along the ten curves 
corresponds to a canard, formed by transversal intersection of the the repelling and attracting 
branches of slow manifolds from figure [7l The curves give the maximum value of x reached along 
each canard solution (compare for example to theoretical results in 28|). The strong canard is 


labeled 7 ^*, and nearby, we see the development of secondary canards, labelled 7 ^ to 7 ®, as k 
increases from zero. 

The number of canards grows very quickly with k at smaller values (A: < 1000). Each new branch 
of secondary canards emerges from one particular curve, labeled as the branch of weak canards, 
7 «)fc branch could not be computed for all values of k, and is partly derived from the envelope 

of endpoints of the secondary canards, however its identification as the weak canard is supported 
for several reasons. Primarily, only the weak canard should coexist with 7 ®* for all parameters, 
and furthermore the bifurcation of secondary canard branches fits with previous results for canards 
in the case of a folded node; we refer the reader to 

expected dia...., and to Q a oo.ptded Mtacatiop dia^tapt. 


28l | for theoretical results and a sketch of the 
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k=im yt=1000 yt=10000 



Figure 7: Slow manifolds in a regularization of the piecewise-smooth system (I20II . plotted in the section {y = 0} for 
different values of the stiffness parameter k. (i) Computations of the slow manifolds, and (ii) magnifications showing the 
increasing numbers of intersections with increasing k. Each intersection of the slow manifolds corresponds to a canard. 
(Colour online indicates the repelling (thicker/blue) and attracting (thinner/red) slow manifolds). 




Figure 8: Branches of canards in a regularization of the piecewise-smooth system (1201) . For valid canard solutions, the 
maximum value of x reached along the canard is plotted against the stiffness parameter k. The strong canard ( 7 “*), weak 
canard ( 7 “*^), and several secondary canards ( 7 ^- 7 ®) are identified. 

Although only the first nine secondary canards are identified here, these computations suggest a 
trend towards infinitely many canards in the nonsmooth (k —?• 00 ) limit. This scenario fits with the 
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model (|26p obtained by pinching, suggesting that the piecewise-smooth system (I26p approaches the 
smooth system m, in the limit where small e leads to an increasingly sharp jump in dynamics near 
the critical manifold. Even allowing for canards beyond those counted here, figure|8]is expected to be 
incomplete, since for large k {k > 1000) additional canard branches were detected in computations, 
born through folds instead of from the branch of weak canards, whose identity it currently unclear. 
One can ask how these various features depend on the form of the sigmoid function that replaces 
sign (17), and whether they are represented in any way in the piecewise-smooth system. These 
are interesting problems for further study, with particular relevence to the study of uniqueness in 
regularization of piecewise-smooth systems. For the present paper, we now return to developing 
the pinch approximation. 


IV. SECOND APPROXIMATION: ROTATING CANARDS 


The exponential scaling in the previous section attempts to resolve the slow dynamics in the 
neighbourhood of the critical manifold u = 0. However, after the microscope on u = 0, the flow 
is still seen to evolve fast towards another surface, the nullcline u = 0. An improvement on this 
method is therefore to take a new microscope and pinch, centred on this new surface. 

The nullcline where ft = 0 is given by Vy in equation (I16j) . Since this is not well defined at z = 0, 
we approximate it by Vq given by equation (I17h . and take a microscope on Vq, by introducing a 
new variable 


u 1 + u 
w = \ --— 

e 2e 


k] 


in terms of which the dynamical system (fT4)) becomes 

w = w {2z + , 

y = 1, 


(28) 


(29) 


i = et^Vk] + . 

This is simulated in figure[9](i). The weak and strong canards now satisfy {tc, z'\ = { —l/(2e)^, /ry/2} 
and {w^z} = {—(^/2e)^, y/2}, respectively. 

The nullcline w = 0 is the curve 


= I (rc, z) G : w = - ' ' | . (30) 

The flow projection onto w = 0 has yet again the same y-z system as before, namely equation (j22jl . 
However, the fact that the flow organises around the nullcline Qy as apposed to Vy gives different 
dynamics, as will be revealed by the pinch. 
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Slow dynamics in figure [9ji) is seen numerically to dominate in an neighbourhood of the 
nullcline h = 0, which corresponds to |ri — < e or |r(;| < 1. In the first approximation of 

section EH slow dynamics is observed to dominate in an e-neighbourhood of the critical manifold, 
\h\ < e, which corresponds to |n| < 1 or |u| < 1. An analytic explanation of why the slow 
neighbourhood is of order e around the critical manifold, and e^ around its associated nullcline, is 
outside the scope of the current paper, but deserves attention in future work. 

The pinch is now enacted similarly to the previous section, by introducing a new variable W = 
w — sign w, which gives 

W = (lT + signVb){2z + ^(lT + signiy)-i/[^l} , 

y = 1, (31) 

i = s{W + signiy)V[d + ^ 

as shown in figure [Hii) . 



Figure 9: The second microscope and pinch, (i) Flow in the microscope system (|29|1 simulated in the plane y = —1 with 
e = 0.05 and yt = 1/8.5. (ii) Flow in the system (I31II obtained by pinching together the surfaces w = ±1 in (i), with the 
vertical axis rescaled for clarity. The primary canards ^Bt,wk jn(jjc-ated. 

The sliding flow is superficially given by the usual equation (I22p . as is found by solving for 
VT = 0 on IT = 0. Crucially, however, we must consider the arrangement of the tangencies W = 0 
on IT = 0, which give the boundaries of the sliding regions. The IT > 0 and IT < 0 subsystems in 
equation (I3ip are tangent to the switching manifold where IT reaches zero as IT approaches zero 
from above or below. These tangencies lie along f = — i^signlT, hence the sliding regions are found 
to be given by 


zjyl > on IT = 0 


(32) 
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illustrated in figure [TO] for different values of jj.. 



*y -►y 


Figure 10: Sliding flow after the second pinch, given by equation (1221) in the regions \zly\ > on IK = 0, for 

£ = 0.05. In (i) there are no canards for ^ = 1/4, and in (ii) the strong canard can be seen in the sliding flow for 
jj, = 1/16. Inset: the strong canard in the microscope system falls inside the pinch region in (ii) but not in (i). 

The curvature of the flow is specified by the second derivative, IT = (1 ± 2e)(/i ± 2e)/2e. Since 
we have 0 < // < 1 and 0 < e <C 1, the flow in IT >0 satisfies IT = (1 + 2e)(/i + 2e)/2e > 0 on its 
tangency line zjy = /i/4e, and hence curves away from the switching manifold. The flow in IT < 0 
satisfies IT = (1 — 2e)(^ — 2e)/2e on its tangency line zjy = —jjLjAe, and hence curves away from 
the switching manifold if < 2e, and towards it if > 2e. Although we are interested in arbitrarily 
small e, either of these can be satished for small enough [i. 

This happens because the weak eigendirection of equation (|22l) lies outside the sliding region. 
The strong eigendirection lies inside the sliding regions if /r > 2e, or outside the sliding regions if 
H < 2 e, meaning that in the former case the sliding dynamics captures no canards, and in the latter 
captures a single canard. (This is immediately in contrast to the continuum of sliding canards in 
system ()27|) ). The two different cases are shown in figure flOl 

The weak and strong canards in the unpinched system lay at {w,z} = { —l/(2e)^,/iy/2} and 
{w,z} = { — {fj,/2£Y,y/2} respectively. The weak canard clearly avoids the pinch region |r(;| < 1 for 
e < 1/2, and then lies at {W,z} = {1 — l/(2e)^,/ry/2}. The strong canard also avoids the pinch 
region if /x > 2e, and is given by {W,z} = {1 — l/(2e)'', y/2}. If /x < 2e the strong canard falls 
inside the pinch region and is not part of system (|31l) for IT 7 ^ 0; if it exists it is part of the sliding 
dynamics on IT = 0. Indeed we see that is exactly the case in hgure fTOl 
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A. Linearizing about the weak canard 

We complete this study by showing that the pinched approximation of the folded node possess 
secondary canards, where (n) denotes the largest integer less than n. These canards rotate 
around the weak canard near y = 0, with rotation numbers taking all integers from 1 to 
then connect to sliding solutions in the attracting and repelling sliding regions that take them to 
y ±oo, as sketched in figure fTTl Note that the rotation takes place in IT <0, therefore secondary 
canards satisfy W < 0. 



Figure 11: Primary canards (dashed) and a secondary canard (bold) with rotation number 3 in the pinched {W,y,z) 
system. 


To solve system (j^TIl we treat the IT > 0 and IT < 0 systems separately, making different 
approximations in the two regions about the dominant singularities. In the region IT < 0, we 
linearize about the weak canard at {W, z} = {1 — {2£)~^, yy/2}. In the region IT > 0, we expand 
about the tangency to the switching manifold at (W, z) = (0, —fiy/Ae). To leading order these give 


in IT > 0 : < 


in IT < 0 : < 


IT 

y 

z 

IT 

y 


M^ + 0(1T), 

1 , 

e + + 0 (IT), 

^ + f (IT - 1 + ^) + 0 {6W^e6z6W) , 

1 , 


(33) 


(34) 


_ i = | + (2e)--i(lT-l + ^) + 0(51T2), 

where 6W = IT — 1 + (l/2e)^ and 6z = z — yy/2. In the remainder of this section we omit the 
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error terms, and find solutions to the truncated local equations. Note that l/(2e)'' ~ 1 —elog(2e) + 
0 ((elog 2 e)^) deviates quickly from unity as e increases from zero, so we cannot approximate it 
by unity. Approximating around the weak canard in Vk < 0 leads to a slight shift in the sliding 
region ()32p . The boundary where the IT > 0 system is tangent to W = 0 is given, as before, by 
zjy = —yLjAe. The boundary where the W < 0 system is tangent to kk = 0 is now given by 

zjy = //(£-! + (2e)^)/2e . (35) 



Figure 12: The pinched system (i) and the expansion p3|) -l|34| ) . Flow simulated in the plane y = —1 for e = 0.05, 
y = 1/8.5. The dotted curves show the nullcline Ik = 0. 


Let US now find the canards in this approximation. We will consider only those that satisfy the 
following properties: 

Cl. y{Q) = 2 ;( 0 ) = 0 . Because canards must satisfy kk —)• 0 as t —f oo±, and the system (I51)) - (I55]1 
is symmetric under the substitution {y, i—>■ {—y, —2, —t}, canards are expected to inherit 

this symmetry, implying that y = z = 0a,tt = 0. 

C2. W{tc) < kk(tc) = kk(tc) = 0 at some t = tc ^ 0- This is because the conditions kk = kk = 0 
define the boundary of the sliding region (equation (|35P ). at which solutions can pass between 
kk < 0 and the invariant sliding region on kk = 0. They can only do so if the solution is 
curving into (kk < 0 ) the region kk < 0 . 


We must then find solutions of the kk < 0 system that satisfy these two conditions. Considering 
system (IM)) we have, noting z = dz/dy = dz/dt, 


1 


z = 


(2s) 


1—e 


»•=- £) 
2e e ' 2' 
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Letting z = (^ + l2 and y = t-\/ ej this rearranges to the Hermite equation 


whose general solution can be written 
C(t) = tC'(O) iFi 1 


C"-rC' + ic = o, 


y. — 1 3 / 1 1 

-> n) TT ) + C(0) 1^1 ( “7^) n) TT ; ) 


(36) 


2 /i ’ 2 ’ 2 


2y T 2 


with derivative 


C(t) = C'(0) iF, ( 


2y ' 2' 2 J y \ 2y 2' 2 


in terms of the confluent hypergeometric function iFi (also known as Rummer’s function M where 
1^1 (a,/3; 7) = M{a,l3,-f), see [l|). 

Applying condition Cl above, at a point where y(0) = 2:(0) = 0 we have VR(0) = 0, hence 
z{0) = 0. In the transformed coordinates this gives initial conditions 

C(o) = 0 , c"(o) = 0 , 


the former of which simplifies the solution above to 


C('^) ^ - 1 3 r" 

= T 


C'W ^ 1 


C'(0) V ’ 2’ 2 y ’ C'(0) " V ' 2’ 2 

or in terms of the Gamma function L and Hermite polynomials Hn [l|, 


(37) 


C('r) = 


*C'(0)e- /"r(i|^) 


il,y,iT) - D+ (0,/i,r)j , 


where 


[m,y,T] = . 

Substituting back in r = tyjy/e^ we find that in LL < 0 there exist solutions given by 

W{t) = l-(2e)--(l + e/r-2£i(0)iFi(^^,l,fg)), 

y{t) = t , 


(38) 


z{t) = ^ + ti(0) iFi 


At —1 3 
2ai ’ 2 ’ 2e 


and these form the portions of any secondary canards that lie in IR <0, outside the switching 
manifold. To this we must apply the second condition, C2, to pick out solutions in system (l38l) that 
tangentially touch (or graze) the boundaries of the sliding regions on VR = 0, where they connect to 
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sliding solutions that form the tails of the canards. Substituting the conditions W(tc) = W(tc) = 0 
at some t = tc ^ 0 into equation (I34p . we find that z{tc) and i(tc) are given by 


z{tc)/[itc = z{tc) + e(l - /r)/2e = (! + £- (2e)^)/2e . (39) 

subject to the boundary conditions 


Figure [T3l shows a simulation of solutions given by system 

(IMD. 

The number and geometry of the secondary canards, specifically the number of rotations they 
make around the weak canard, are easily found as follows. The confluent hypergeometric function 
iFi{—a,b,c) has 2{a + 1) real zeros [l| (with (n) denoting the largest integer smaller than n), 
between which the function oscillates through 2(a + 1) — 1 maxima/minima, and between these the 
functions makes (a) complete oscillations. These oscillations form the rotations of the secondary 
canards. 

Using the boundary conditions W{tc) = W{tc) = 0, a given value of tc > 0 picks out one 
of the solutions (I38p . with rotations. The boundary conditions may be satished at any 

one of the maxima, so any of the rotation numbers from 1 up to are obtained, by different 


solutions with unique values of tc- Hence there exist {-^) secondary canards with rotation numbers 



Figure 13: Secondary canards for e = 0.05 and ^ = 1/8.5, from the Hermite solution (13811 . Inset: the Hermite solution 
(full curves) is compared to the approximation using (|40p (dotted curves). Long dashes indicate segments of sliding. The 
weak and strong canards are shown. 


As a final remark, the complicated expressions above can be simplified with the large parameter 
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asymptotic approximation 

iFi {abx) = ^ (f (6 - cos (^^2x{b-2a) + f (1 - 26)) (40) 

(see e.g. Eq.13.5.14 of [l|) This is sufficient to approximate solutions for small x = but this 

exponential approximation (compared to the exact solution in figure fTSDiil) is not accurate enough 
to correctly give the correct number of, or number of rotations of, secondary canards. 


V. CLOSING REMARKS 


Pinching captures the key geometry - singularities and bifurcations - necessary to provide a 
discontinuous model of a singularly perturbed dynamical system. As a method of approximation it 
is purely qualitative, yet it appears to accurately describe singular features such as both primary 
and secondary canards in the folded node studied here. In a previous study of the van der Pol 
oscillator Q], pinching was also shown to capture the maximal canard, and to distinguish between 
a canard explosion and a Hopf bifurcation. 

Our main aim with this work is to help illuminate the bridge between smooth and piecewise- 
smooth models of dynamical systems, by showing how closely the phenomena of singular perturba¬ 
tion and discontinuity are related. Canards in slow-fast systems have now been known for more than 
three decades, having first been studied using non-standard analysis j4|, and later with standard 


tools such as matched asymptotic expansions and geometric singular perturbation theory Q,Q 
Canards in discontinuous system lay unannounced in the work of Filippov [l^, 171 for considerably 
longer. Only recently has the link between folds in critical manifolds of slow-fast systems, and the 
two-fold singularity in discontinuous systems, become clear, through the methods of regularization 
23 1 and of pinching Q], the former applying a topological equivalence in the singular limit e = 0 , 


the latter approximating the geometry for e 7 ^ 0 . 
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